Inelastic semiclassical Coulomb scattering 



Gerd van de Sandf and Jan M RostJ 

f— Theoretical Quantum Dynamics - 

Fakultat fiir Physik, Universitat Freiburg, Hermann-Herder— Str. 3, D-79104 
Freiburg, Germany 

JMax-Planck-Institute for Physics of Complex Systems, Nothnitzer Str. 38, 
D-01187 Dresden, Germany 



Abstract. We present a semiclassical S-matrix study of inelastic collinear 
electron-hydrogen scattering. A simple way to extract all necessary information 
from the deflection function alone without having to compute the stability matrix 
is described. This includes the determination of the relevant Maslov indices. 
Results of singlet and triplet cross sections for excitation and ionization are 
reported. The different levels of approximation - classical, semiclassical, and 
uniform semiclassical - are compared among each other and to the full quantum 
result. 



PACS numbers: 34.80D, 03.65. Sq, 34.10+x 

1. Introduction 

Semiclassical scattering theory was formulated almost 40 years ago for potential 
scattering in terms of WKB-phaseshifts Ten years later, a multidimensional 

formulation appeared, derived from the Feynman path integral j^j . Based on a similar 
derivation Miller developed at about the same time his 'classical S-matrix' which 
extended Pechukas' multidimensional semiclassical S-matrix for potential scattering 
to inelastic scattering || |J, |j| . These semiclassical concepts have been mostly applied 
to molecular problems, and in a parallel development by Balian and Bloch Q to 
condensed matter problems, i.e. to short range interactions. 

Only recently, scattering involving long range (Coulomb) forces has been 
studied using semiclassical S-matrix techniques, in particular potential scattering 0, 
ionization of atoms near the threshold || and chaotic scattering below the ionization 
threshold pc|| . The latter problem has also been studied purely classically jll] and 
semiclassically within a periodic orbit approach |T^| . 

While there is a substantial body of work on classical collisions with Coulomb 
forces using the Classical Trajectory Monte Carlo Method (CTMC) almost no 
semiclassical studies exist. This fact together with the remarkable success of 
CTMC methods have motivated our semiclassical investigation of inelastic Coulomb 
scattering. To carry out an explorative study in the full (12) dimensional phase space 
of three interacting particles is prohibitively expensive. Instead, we restrict ourselves 
to collinear scattering, i.e. all three particles are located on a line with the nucleus 
in between the two electrons. This collision configuration was proven to contain the 
essential physics for ionization near the threshold [§, [l3[ Q and it: nts wel1 into tnc 



context of classical mechanics since the collinear phase space is the consequence of a 
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stable partial fixed point at the interelectronic angle #12 = 180° Moreover, it is 
exactly the setting of Miller's approach for molecular reactive scattering. 

For the theoretical development of scattering concepts another Hamiltonian of 
only two degrees of freedom has been established in the literature, the s-wave model 
p5[ . Formally, this model Hamiltonian is obtained by averaging the angular degrees 
of freedom and retaining only the zeroth order of the respective multipole expansions. 
The resulting electron-electron interaction is limited to the line n = ri , where the n 
are the electron-nucleus distances, and the potential is not differentiable along the line 
r\ =T2- This is not very important for the quantum mechanical treatment, however, 
it affects the classical mechanics drastically. Indeed, it has been found that the s-wave 
Hamiltonian leads to a threshold law for ionization very different from the one resulting 
from the collinear and the full Hamiltonian (which both lead to the same threshold 



law) 16 1. Since it is desirable for a comparison of semiclassical with quantum results 
that the underlying classical mechanics does not lead to qualitative different physics 
we have chosen to work with the collinear Hamiltonian. For this collisional system we 
will obtain and compare the classical, the quantum and the primitive and uniformized 
semiclassical result. For the semiclassical calculations the collinear Hamiltonian was 
amended by the so called Langer correction, introduced by Langer ]l7f to overcome 
inconsistencies with the semiclassical quantization in spherical (or more generally non- 
cartesian) coordinates. 

As a side product of this study we give a rule how to obtain the correct Maslov 
indices for a two-dimensional collision system directly from the deflection function 
without the stability matrix. This does not only make the semiclassical calculation 
much more transparent it also considerably reduces the numerical effort since one 
can avoid to compute the stability matrix and nevertheless one obtains the full 
semiclassical result. 

The plan of the paper is as follows: in section 2 we introduce the Hamiltonian and 
the basic semiclassical formulation of the S-matrix in terms of classical trajectories. 
We will discuss a typical S-matrix S(E) at fixed total energy E and illustrate a simple 
way to determine the relevant (relative) Maslov phases. In section 3 semiclassical 
excitation and ionization probabilities are compared to quantum results for singlet and 
triplet symmetry. The spin averaged probabilities arc also compared to the classical 
results. In section 4 we will go one step further and uniformize the semiclassical S- 
matrix, the corresponding scattering probabilities will be presented. We conclude the 
paper with section 5 where we try to assess how useful semiclassical scattering theory 
is for Coulomb potentials. 

2. Collinear electron-atom scattering 

2.1. The Hamiltonian and the scattering probability 

The collinear two-electron Hamiltonian with a proton as a nucleus reads (atomic units 
are used throughout the paper) 

h = & + &-L-L--*—. (i) 

2 2 n r2 n + r 2 

The Langer-corrected Hamiltonian reads 

H = h+ ii + k < 2 > 
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For collinear collisions we have only one 'observable' after the collision, namely the 
state with quantum number n, to which the target electron was excited through the 
collision. If its initial quantum number before the collision was ml , we may write the 
probability at total energy E as 

P n , n ,(E) = \(n\S\n')\ 2 (3) 

with the S-matrix 

S= Ihn e «/* e -^(*-*') e - iff t*'. (4) 

t'— * — oo 

Generally, we use the prime to distinguish initial from final state variables. The 
Hamiltonians Hi and Hf represent the scattering system before and after the 
interaction and do not need to be identical (e.g. in the case of a rearrangement 
collision). The initial energy of the projectile electron is given by 

e' = E-e' (5) 

where e' is the energy of the bound electron and E the total energy of the system. 
In the same way the final energy of the free electron is fixed. However, apart from 
excitation, ionization can also occur for E > in which case \n) is simply replaced by 
by a free momentum state \p). This is possible since the complicated asymptotics of 
three free charged particles in the continuum is contained in the S-matrix. 



2.2. The semiclassical expression for the S-matrix 
Semiclassically, the S-matrix may be expressed as 



V'(s)=Ef2.(s)^ rif ^ (6) 

3 

where the sum is over all classical trajectories j which connect the initial state n' and 
the final 'state' n with a respective probability of V^ n ,(E). The classical probability 

Pnl(E) ^ given by 



n,n'\ ) e , e ' V ) Q n N 



de(R') 



see H where also an expression for the normalization constant N is given. Note, that 
due to the relation (||) derivatives of e and e with respect to n or R' differ only by a 
sign. From now on we denote the coordinates of the initially free electron by capital 
letters and those of the initially bound electron by small letters. If the projectile is 
bound after the collision we will call this an 'exchange process', otherwise we speak of 
'excitation' (the initially bound electron remains bound) or ionization (both electrons 
have positive energies). The deflection function e(R') has to be calculated numerically, 
as described in the next section. The phase <frj is the collisional action given by 

$j- (P, n; P', n') = - J dt (qn + RP) (8) 

with the angle variable q. The Maslov index Vj counts the number of caustics 
along each trajectory. 'State' refers in the present context to integrable motion for 
asymptotic times t — > ±oo, characterized by constant actions, J' = 2i:h(n' + 1/2). 
The (free) projectile represents trivially integrable motion and can be characterized by 
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Figure 1. Scattering trajectories at a total energy of E = 0.125 a.u. with initial 
conditions marked in figure El The lahcls (a-f) refer to representative trajectories 
with initial values R' shown in figure bl The left column corresponds to classical 
exchange n' = 1 — > n = 1, the middle column represents ionization events and 
the right column shows elastic back-scattering with n' = 1 — * n = 1. 



its momentum P' . In our case, each particle has only one degree of freedom. Hence, 
instead of the action J' we may use the energy e' for a unique determination of the 
initial bound state. In the next sections we describe how we calculated the deflection 
function, the collisional action and the Maslov index. 



2.2.1. Scattering trajectories and the deflection function The crucial object for the 
determination of (semi-)classical scattering probabilities is the deflection function 
e(R') where e is the final energy of the projectile electron as a function of its initial 
position Rq + R' . Each trajectory is started with the bound electron at an arbitrary 
but fixed phase space point on the degenerate Kepler ellipse with energy e' = —1/2 
a.u.. The initial position of the projectile electron is changed according to R' , but 
always at asymptotic distances (we take Rq = 1000 a.u.), and its momentum is fixed 
by energy conservation to P 1 = [2(E — e')] 1 / 2 . The trajectories are propagated as a 
function of time with a symplectic integrator []l9f and e = e(t — > oo) is in practice 
evaluated at a time t when 

dln\e\/dt<5 (9) 

where S determines the desired accuracy of the result. Typical trajectories are shown 
in figure [I], their initial conditions are marked in the deflection function of figure |^. 

In the present (and generic) case of a two-body potential that is bounded from 
below the deflection function must have maxima and minima according to the largest 
and smallest energy exchange possible limited by the minimum of the two-body 
potential. The deflection function can only be monotonic if the two-body potential is 
unbounded from below as in the case of the pure (homogeneous) Coulomb potential 
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Figure 2. The deflection function at an energy of E = 0.125 a.u. and for an 
initial state as described in the text. The energy interval enclosed by dashed lines 
marks ionizing initial conditions and separates the exchange region (e < 0) from 
the excitation region (e > E), where e is the energy of the projectile after the 
collision. 



without Langer correction (compare, e.g., figure 1 of ||). This qualitative difference 
implies another important consequence: For higher total energies E the deflection 
function is pushed upwards. Although energetically allowed, for E > 1 a.u. the 
exchange-branch vanishes as can be seen from figure ||. As we will see later this 
has a significant effect on semiclassical excitation and ionization probabilities. 

2.2.2. The form of the collisional action The collisional action <!>., along the 
trajectory j in (j^) has some special properties which result from the form of the 
S-matrix (|4|). The asymptotically constant states are represented by a constant action 
J or quantum number n and a constant momentum P for bound and free degrees of 
freedom respectively. Hence, in the asymptotic integrable situation with h = P = 
before and after the collision no action is accumulated and the collisional action has 
a well defined value irrespectively of the actual propagation time in the asymptotic 
regions. This is evident from (|J) which is, however, not suitable for a numerical 
realization of the collision. The scattering process is much easier followed in coordinate 
space, and more specifically for our collinear case, in radial coordinates. In the 
following, we will describe how to extract the action according to (||) from such a 
calculation in radial coordinates (position r and momentum p for the target electron, 
R and P for the projectile electron). The discussion refers to excitation processes 
to keep the notation simple but the result (O) holds also for the other cases. The 
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Figure 3. The deflection function at an energy of E = 2 a.u. and for an initial 
state as described in the text. The dashed line separates ionizing initial conditions 
from excitation events. 



collisional action <!> can be expressed through the action in coordinate space $ by Q 
$(P,n;P',n') = $(P, r; P' , r') + F 2 (r' , ri) - F 2 (r, n), (10) 

where 



r; P' 



= lim 



dr 



-RP + pr 



(11) 



is the action in coordinate space and F 2 is the generator for the classical canonical 
transformation from the phase space variables (r,p) to (q,n) given by 



F 2 (r, n) — sgn(p) / (2m [e (n) — v (x)]) 2 dx . 



(12) 



Here, denotes an inner turning point of an electron with energy e(n) in the potential 
v(x). Clearly, F 2 will not contribute if the trajectory starts end ends at a turning point 
of the bound electron. Partial integration of ( |TT| ) transforms to momentum space and 
yields a simple expression for the collisional action in terms of spatial coordinates: 



$(P,n;P',n') = lim -J dr RP 



rp 



(13) 



Note, that t[ and refer to times where the bound electron is at an inner turning 
point and the generator F 2 vanishes. Phases determined according to ( [j"3| ) may still 
differ for the same path depending on its time of termination. However, the difference 
can only amount to integer multiples of the (quantized !) action 

J = <j> pdr = 2n\n + )\ (14) 

of the bound electron with e < 0. Multiples of 2-k for each revolution do not change 
the value of the S-matrix and the factor ^ is compensated by the Maslov index. In 
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the case of an ionizing trajectory the action must be corrected for the logarithmic 
phase accumulated in Coulomb potentials |lS{ |. 

Summarizing this analysis, we fix the (in principle arbitrary) starting point of 
the trajectory to be an inner turning point (r'^p' = 0,p' > 0) which completes the 
initial condition for the propagation of trajectories described in section 2.2.1. In 
order to obtain the correct collisional action (||) in the form (|l3|) we also terminate a 
trajectory at an inner turning point r% after the collision such that $ is a continuous 
function of the initial position R' . Although this is not necessary for the primitive 
semiclassical scattering probability which is only sensitive to phase differences up 
to multiples of J as mentioned above, the absolute phase difference is needed for a 
uniformized semiclassical expression to be discussed later. 




2.3. Maslov indices 

2.3.1. Numerical procedure In position space the determination of the Maslov index 
is rather simple for an ordinary Hamiltonian with kinetic energy as in (|^). According 
to Morse's theorem the Maslov index is equal to the number of conjugate points along 
the trajectory. A conjugate point in coordinate space is defined by (/ degrees of 
freedom, (qi,Pi) a pair of conjugate variables) 

det {M qp ) = del I I = . (15) 

^ (//,.- ' ^ 

The matrix M qp is the upper right part of the stability or monodromy matrix which 
is defined by 

In general, the Maslov index Uj in @ must be computed in the same representation 
as the action. In our case this is the momentum representation in (|l3|). However, 
the Maslov index in momentum space is not simply the number of conjugate points 
in momentum space where det (M pq ) — 0. Morse's theorem relies on the fact that 
in position space the mass tensor By = d 2 H/dpidpj is positive definite. This is not 
necessarily true for Dij = d 2 H/dqidqj which is the equivalent of the mass tensor in 
momentum space. How to obtain the correct Maslov index from the number of zeros of 
det (M pq ) = is described in |pp[ , a review about the Maslov index and its geometrical 
interpretation is given in pl[ . 



2.3.2. Phenomenological approach for two degrees of freedom For two degrees of 
freedom, one can extract the scattering probability directly from the deflection 
function without having to compute the stability matrix and its determinant explicitly 
In view of this simplification it would be desirable to determine the Maslov 
indices also directly from the deflection function avoiding the complicated procedure 
described in the previous section. This is indeed possible since one needs only the 
correct difference of Maslov indices for a semiclassical scattering amplitude. 

A little thought reveals that trajectories starting from branches in the deflection 
function of figure separated by an extremum differ by one conjugate point. This 
implies that their respective Maslov indices differ by Av = 1. For this reason it is 
convenient to divide the deflection function in different branches, separated by an 
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extremum. Trajectories of one branch have the same Maslov index. Since there are 
two extrema we need only two Maslov indices, vi = 1 and v 2 = 2. The relevance 
of just two values of Maslov indices (1,2) can be traced to the fact that almost all 
conjugate points are trivial in the sense that they belong to turning points of bound 
two-body motion. 

We can assign the larger index v 2 — 2 to the trajectories which have passed one 
more conjugate point than the others. As it is almost evident from their topology, 
these are the trajectories with de/dR' > shown in the upper row of figure |l|. (They 
also have a larger collisional action <&j). The two non-trivial conjugate points for these 
trajectories compared to the single conjugate point for orbits with initial conditions 
corresponding to de/dR' < can be understood looking at the ionizing trajectories (b) 
and (e) of each branch in figure |]. Trajectories from both branches have in common 
the turning point for the projectile electron (P = 0). For trajectories of the lower row 
all other turning points belong to complete two-body revolutions of a bound electron 
and may be regarded as trivial conjugate points. However, for the trajectories from 
the upper row there is one additional turning point (see, e.g., figure |l|(b)) which cannot 
be absorbed by a complete two-body revolution. It is the source for the additional 
Maslov phase. 

We finally remark that de/dR' > is equivalent to dn/dq < of |^5| leading to 
the same result as our considerations illustrated above. 



3. Semiclassical scattering probabilities 

Taking into account the Pauli principle for the indistinguishable electrons leads to 
different excitation probabilities for singlet and triplet, 

P+(E) = \S e ^(E) + S E -^(E)\ 2 

P-(E)= \S £ , £ ,(E)-S E - e , e ,(E)\ 2 , (17) 

where the probabilities are symmetrized a posteriori (see ]2^| ). Here, S et e' denotes 
the S-matrix for the excitation branch, calculated according to while SE-e,e' 
represents the exchange processes, at a fixed energy e < 0, respectively. 

Ionization probabilities are obtained by integrating the differential probabilities 
over the relevant energy range which is due to the symmetrization (^|) reduced to 
E/2: 

E/2 

P?on( E )= J Pt{E)de. (18) 
o 



3.1. Ionization and excitation for singlet and triplet symmetry 

We begin with the ionization probabilities since they illustrate most clearly the 
effect of the vanishing exchange branch for higher energies as illustrated in figure ||. 
The semiclassical result for the Langer Hamiltonian (^) shows the effect of the 
vanishing exchange branch in the deflection function figure ^ which leads to merging 
P^ probabilities at a finite energy E, in clear discrepancy to the quantum result, 
see figure [|. Moreover, the extrema in the deflection function lead to the sharp 
structures below E = 1 a.u.. The same is true for the excitation probabilities where 
a discontinuity appears below E = 1 a.u. (figure p3). Note also that due to the 
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Figure 4. Ionization probabilities for singlet and triplet according to ( |18[ ) with 
the Hamiltonian (0) (solid line) compared to quantum mechanical calculations 
(dotted line). 
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Figure 5. Semiclassical excitation probabilities (n = 1, 2, 3) according to ( Jl 7J ) for 
singlet (part a) and triplet (part b) in the LSA (solid line) compared to quantum 
mechanical calculations (dotted line). 



violated unitarity in the semiclassical approximation probabilities can become larger 
than unity, as it is the case for the n = 1 channel. 

Singlet and triplet excitation probabilities represent the most differential 
scattering information for the present collisional system. Hence, the strongest 
deviations of the semiclassical results from the quantum values can be expected. Most 
experiments do not resolve the spin states and measure a spin-averaged signal. In our 
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Figure 6. Spin averaged quantum results for ionization (dotted line) compared 
to averaged semiclassical probabilities (solid line) from ( |19[ ) and classical 
probabilities (dashed line) from (E0|). 



model this can be simulated by averaging the singlet and triplet probabilities to 

P € {E)= 1 -{P+{E)+P-{E)). (19) 

The averaged semiclassical probabilities may also be compared to the classical result 
which is simply given by 

P? L {E) = J2CPB(E)+V { S_ e ,(E)) (20) 

3 

with V { Jl{E) from @). 

Figure || shows averaged ionization probabilities. They are very similar to each 
other, and indeed, the classical result is not much worse than the semiclassical result. 

In figure ^ we present the averaged excitation probabilities. Again, on can see the 
discontinuity resulting from the extrema in the deflection function. As for ionization, 
the spin averaged semiclassical probabilities (figure ^d) are rather similar to the 
classical ones (figure ^), in particular the discontinuity is of the same magnitude 
as in the classical case and considerably more localized in energy than in the non- 
averaged quantities of figure ||. 

Clearly, the discontinuities are an artefact of the semiclassical approximation. 
More precisely, they are a result of the finite depth of the two-body potential in the 
Langer corrected Hamiltonian (^|). Around the extrema of the deflection function 
the condition of isolated stationary points, necessary to apply the stationary phase 
approximation which leads to (||), is not fulfilled. Rather, one has to formulate a 
uniform approximation which can handle the coalescence of two stationary phase 
points. 

4. Uniformized scattering probabilities 

We follow an approach by Chester et. al. [^3) . The explicit expression for the uniform 
S-matrix goes back to Connor and Marcus p3] who obtained for two coalescing 
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Figure 7. Spin averaged quantum results (dotted line) for excitation (n = 1, 2, 3) 
compared to classical probabilities (solid line, part a) from (^) and averaged 
semiclassical probabilities (solid line, part b) from (fly). 



trajectories 1 and 2 



S n , n ,(E) = Bi+ {-z) ^U^e^+^+Br (-z) yJv™,(E) 

where 

Bi* (-z) = V¥ \z$Ai {-z) T iz~* Ai' (-*)] e ±l ( 3Z2 ~~) (22) 

2 

The argument z = [| ($2 — ^1)] 3 of the Airy function Ai(z) contains the absolute 
phase difference. We assume that $2 > $1 which implies for the difference of the 
Maslov indices that ^2 — v\ = 1 (compare (^ with (]2l]) and (p3[)). Since the absolute 
phase difference enters ( ^l| ) it is important to ensure that the action is a continuous 
function of R! avoiding jumps of multiples of 2ir, as already mentioned in section 2.2.2. 
For large phase differences (||) is recovered since 

lim Bi ± (-z) = 1 . (23) 



Our uniformized S-matrix has been calculated by applying (|21]) to the two 
branches for exchange and excitation separately and adding or subtracting the results 
according to a singlet or triplet probability. In the corresponding probabilities 
of figure H the discontinuities of the non-uniform results are indeed smoothed in 
comparison with figure ||. However, the overall agreement with the quantum 
probabilities is worse than in the non-uniform approximation. A possible explanation 
could lie in the construction of the uniform approximation. It works with an integral 
representation of the S-matrix, where the oscillating phase (the action) is mapped 
onto a cubic polynomial. As a result, the uniformization works best, if the deflection 
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Figure 8. Uliiformized semiclassical excitation probabilities (n = 1,2,3) 
according to (Ell) (solid line) for singlet (part a) and triplet (part b) compared to 
quantum mechanical calculations (dotted line). 



function can be described as a quadratic function around the extremum. Looking at 
figure ^ one sees that this is true only in a very small neighborhood of the extrema 
because the deflection function is strongly asymmetric around these points. We also 
applied a uniform approximation derived by Miller which gave almost identical 
results. 

Finally, for the sake of completeness, the spin averaged uniform probabilities 
are shown in figure ^. As can be seen, the discontinuities have vanished almost 
completely. However, the general agreement with quantum mechanics is worse than for 
the standard semiclassical calculations, similarly as for the symmetrized probabilities. 

5. Conclusion 

In this paper we have described inelastic Coulomb scattering with a semiclassical S- 
matrix. To handle the problem for this explorative study we have restricted the phase 
space to the collinear arrangement of the two electrons reducing the degrees of freedom 
to one radial coordinate for each electron. In appreciation of the spherical geometry we 
have applied the so called Langer correction to obtain the correct angular momentum 
quantization. Thereby, a lower bound to the two-body potential is introduced which 
generates a generic situation for bound state dynamics since the (singular) Coulomb 
potential is replaced by a potential bounded from below. The finite depth of the two- 
body potential leads to singularities in the semiclassical scattering matrix (rainbow 
effect) which call for a uniformization. 

Hence, we have carried out and compared among each other classical (where 
applicable), semiclassical, and uniformized semiclassical calculations for the singlet, 
triplet and spin-averaged ionization and excitation probabilities. Two general trends 
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Figure 9. Spin averaged uniformizcd excitation probabilities (n = 1, 2, 3, solid 
line) compared to quantum results (dotted line). 



may be summarized: Firstly, the simple semiclassical probabilities are overall in 
better agreement with the quantum results for the singlet /triplet observables than the 
uniformized results. The latter are only superior close to the singularities. Secondly 
for the (experimentally most relevant) spin-averaged probabilities the classical (non- 
symmetrizable) result is almost as good as the semiclassical one compared to the 
exact quantum probability. This holds for excitation as well as for ionization. Hence, 
we conclude from our explorative study that a full semiclassical treatment for spin- 
averaged observables is probably not worthwhile since it does not produce better 
results than the much simpler classical approach. Clearly, this conclusion has to be 
taken with some caution since we have only explored a collincar, low dimensional 
phase space. 
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